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A possible SU(15) process at HERA is investigated. The process that we consider is 
e~P — > u e /i~ fi~ + anything through the exchange of new heavy gauge bosons X~~ and 
X which are predicted in SU(15). This process produces two easily observable like- 
sign muons in the final state. The cross section of this process is calculated by using 
HELAS and VEGAS programs, and PDF-library functions. The cross section turns out 
to be small to be observed in near future. 



The standard model, SU(3)<g>SU(2)<g)U(l) for strong and electroweak interac- 
tions, has been remarkably successful in particle physics. In this model, however, 
there are too many coupling constants. A straightforward extension of this model 
is SU(5). The critical drawback of this model is, however, the nonobservation of 
the proton decay predicted i». the model. There have been many grand unifica- 
tion models including SU(15).EJ We have no profound reason to prefer one model to 
others. 

SU(15) is a nonrenormalizable theory since there is no anomaly cancellation in 
itself. One way to avoid the nonrenormalizability is to introduce mirsor fermion. 
Such introduction of mirror fermion is not favorable but not excluded!!! Therefore, 
it is worth reviewing the possible processes arising from SU(15). 

In this letter, we investigate the possibility to discover the like-sign dilepton 
production, e~ P — > v e l~l~ + anything, in SU(15) in the future e~P collider. The 
dileptons are produced via the exchange of new gauge bosons and they can be 
electrons, muons, or taus. Here, we consider two muons to compare our result to 
that of Agrawal et alu p 

We used the HELAS program!] in the numerical calculation of the cross section of 
the process and the VEGAS programQ for the multi-dimensional (10-dimensional in 
this case) integrations stemmed from the phase spaces of the kinematic variables of 
the particles, v e , pT 's and the EHLQ2 nucleon structure function setu for the d quark 
structure function. The HELAS program is designed mainly to provide numerical 
means for the calculation of cross sections in the standard model. However, it is 
easy to implement new interactions from extended standard models like SU(15). 

The process that we considered is e~ P — > v e u~ u~ + anything through the ex- 
change of heavy gauge bosons X~ or A . We make use of the dilepton interaction 
Lagrangian given in the paper of Frampton et aln 
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+ f|^^(^)-e + ^-^(^)^. (1) 

Where, X~(X + ) and X~~ (X ++ ) are the massive lepton gauge bosons which carry 
two lepton numbers and result from the SU(3)/ breaking. The coupling constant 
33/ is given approximately by 1.19e. Notice that vector couplings are missing for 
the currents with the doubly charged gauge bosons due to Fermi statistics. 

For the process, e~ P — > D e j-i~ (i~ + anything, in e~ P collider, we are interested 
in the subprocess, ed — > v e \i~ fJ,~u, through the exchange of gauge bosons, X's. 
Where, d quark in the process comes as a parton inside the proton. This process 
produces two muons and an electron- neutrino, P e , in the final state while the only 
incoming lepton is the electron. Therefore, although the total lepton numbers are 
conserved, lepton numbers in each generation are violated. Since the standard 
model is free from this violation, this process produces a clean signal without such 
background processes in the standard model. The relevant Feynman graphs for the 
process are shown in figs. 1-4. 

The new gauge bosons can be produced from the incoming lepton. Then, they 
produce two muons directly (as shown in fig. 4) or interact with W boson to produce 
two muons (as shown in figs. 2 and 3) or interact with W boson to produce another 
doubly charged gauge boson which then turns into two muons (as shown in fig. 1). 
The corresponding amplitudes are obtained from the Lagrangian given in eq. (|l|), 
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Where, [i\ and [i 2 indicate two muons. The vertex, WX~X , is similar to triple 
gauge boson vertices like WWW or ggg in the standard model. U u d is the Cabbibo- 
Kobayashi-Maskawa (CKM) matrix element for u and d quarks. Here, we have used 
unitary gauge for the gauge bosons to make use of HELAS subroutines which are 
written in unitary gauge when the bosons are massive. Only a small modification 
such as the change in couplings is needed to apply HELAS subroutines to this 
extended standard model. The total amplitude for the process e~ d — > v e \i~ fi~u is 
the sum of all the amplitudes in eq. (0), 



M = M(l) + M(2) + M (3) + M(4). 



(3) 



We can simplify the differential cross section for the process by writing the four 
body phase space as a product of a few two body phase spaces. 



do 
d 4 (PS) 



-di(PS), 



2(2n) s \ 1 / 2 (s,m e ,m d ) 
d 2 (e + d -> A + B)dm 2 A dm 2 B d2{A 



v e u)d 2 {B — > n (j, ). (4) 



Where, da is the differential cross section for a parton (quark) q from the proton 
to be scattered by the electron, s — xs is the square of the center of mass energy of 
the electron and the parton from the proton with the total center of mass energy 
y/s of the system of the incoming electron and the incoming proton, x is the 
standard Bjorken variable and thus the d quark carries the fractional momentum 
of the proton, x times the proton momentum. The four body phase space d±{PS) 
can be obtained from the process e~d — > i> e fi~ fj,~u via the process e _ + d — > A + B 
and then A — > v e + u and B — * /i~ + The two body phase spaces are given as 



d 2 (e + d-+ A + B) 
d 2 {A — * v e u) 
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m,A is the invariant mass of the system of the u quark and the antincutrino and ms 
is the invariant mass of the system of two muons. Angular variables are defined in 
their center of mass frames 

Similarly, c?f2* is defined in the center of mass frame of u and u, 

of two muons, [i~ and fi~ . 

We can calculate the total cross section by summing over all the partons and 
integrating over x. 



dQ* B is defined in the center of mass frame of A and B. 

and dQ*. in that 



j dxy2F q (x,Q 2 )o(V§ = xs,M x ). 
Jo 
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Here, instead of summing over all flavors, we need only d quark distribution function 
to produce two like-sign muons and an antineutrino in the final state. Therefore, the 
summation of the distribution functions over all the flavors can be simply replaced 
by the d quark distribution function, '}2 q F q (x 1 Q 2 ) — » d(x,Q 2 ). We use this d 
quark structure function from the EHLQ2 nucleon structure function set in the 
PDF- library routines in the Phenomenology Institute at the Univesity of Wisconsin. 
Taking different nucleon structure function sets such as EHLQ1 or KMRSB does 
not make any physically significant difference. 

The numerical value of the cross section of the process is obtained by using 
HELAS and VEGAS programs. The VEGAS adaptive monte carlo integration 
subroutines were used for the multi-dimensional integrations arising from the kine- 
matics of the multi-particles involved in the process. In the VEGAS program, the 
'FUNCTION' subprogram for the cross section has to be defined in connection with 
the HELAS program. This FUNCTION subprogram represents the cross section 
and thus the HELAS program plays the central role in FUNCTION. 

Before starting HELAS, we have to prepare a few subroutines. For the coupling 
which involve in the interactions of the process, the relevant coupling subroutines 
(COUP IX and COUP2X) should be determined. To find the momenta of the par- 
ticles in the center of mass frame, we use the momentum subroutine (MOM2CX) 
along with the subroutines (MOMNTX, BOOSTX) related to the boost factors. 
Now, we draw the Feynman diagrams of this process and then apply HELAS Feyn- 
man rules. This is the central part of the application of HELAS. After setting up 
all these, we are ready to use HELAS and follow the steps below. 

First, we find the wave functions (IXXXXX, OXXXXX) of the external particles 
by identifying flowing in and flowing out particles (all are fermions in this case). 
Second, we find the vertex subroutines ( JIOXXX, FVOXXX, VVVWXX, IOVXXX) 
for the corresponding vertices in the Feynman diagrams. Third, we calculate the 
total amplitude with the phase space factors. Fourth, we multiply EHLQ2 u quark 
structure function_(EHLQ2PDF). Agrawal et al. used EHLQ1 for the u quark 
structure function.EI To obtain the numerical result for the cross section, we have to 
integrate over all relevant kinematic variables by running VEGAS with the above 
FUNCTION subprogram. 

In e~ P collider at HERA, electrons of the energy 30GeV collide on protons of the 
energy 820GeV in producing the center of mass energy, y/s = 314GeV. We chose a 
few values of the mass of the gauge boson, M x «100GeV, 200GeV, or 300GeV. We 
also used the constituent quark masses, m u = 0.3GeV and = 0.3GeV. All other 
standard values were used in the numerical calculation such as Mw = 80.22GeV 
and sinOw = 0.2325GeV, etc. With all the values considered, we obtained the 
following results for the cross section (see Table 1). 

The cross section of the process is 8 x l(T 2 fb for M x = fOOGeV, 3 x I CT 4 fb 
for M x = 200GeV, and 3.5 x l(T 6 fb for M x = 300GeV. These results show that 
the cross section decreases very sharply as the gauge boson X increases. Since the 
cross section is very small, it is quite difficult to observe the process in the fore- 
seeable future. The result is also an order of-,2 smaller than the one in the process 
considered by Agrawal, Frampton, and Ng.H There, they considered the process, 
e~q — ► e + 2/i~q. In the process, there is no quark flavor changes and there is a 
positron instead of antineutrino in the final state. The behavior of the sharp dec- 
lination of the cross section as the increment of M x is similar to both our process 
and theirs. Although the cross section of our process is much smaller than that of 
theirs, our process does not contain well measurable particle, e + , in the final state 
with different kinematics. In our process, there is the missing momentum due to 
neutrinos without electron or positron in the final state. Therefore, we hope to be 
able to isolate our process from theirs in the collider of next generation. Unfortu- 
nately, however, considering the current luminocity ~ 10 31 cm _2 s _1 of HERA we 
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expect that the cross section of this process is too small to be observed in near 
future. g 

Frampton and Ng3 showed that the lower bound of the gauge boson mass, Mx , 
is greater than 120GeV from the experimental .angular distribution in Bhabha scat- 
tering, e + e~ — > e + e~. Carlson and FramptonQ concluded Mx >230GeV from the 
precision data of the angular distribution of electrons in polarized muon decay. Since 
the value of the coupling, g%i ss 1.19e is not the unavoidable choice, the above mass 
lower bound is not strict. And thus Agrawal et aln calculated the cross section of 
the process, e~q — ► e + 2fi~q, with the broad range of Mx, lOOGeV < Mx < ITeV. 
We also calculated the cross section of the process, e~ P — > v e /j,~ fjT + anything, 
starting with Mx =100GeV to compare our result to that of Agrawal et aln. How- 
ever, the lower bound of Mx should be well above lOOGeV. Another good way to 
find the mass limit of the dilepton gauge boson may be to look into the precision 
electroweak data at the Z peak. The loop calculation with these dilepton gauge 
bosons should modify the radiative correction to Z or W boson propagators. The 
comparision of the result with the current data may set a better lower limit on 
Mx- We have to keep in mind that SU(15) is not a renormalizable theory and some 
modifications have to be made in the application of the HELAS program in this 
calculation. 

Table 1. The cross section for e~ + d — > i> + u + /x~ + fi~ with various value of Mx 



Mx(GeV) 


cross section in 


(fb) 


cross section of AFNH in (pb) 


100 


a ~ 8.0 x 10" 


-la 


a ~ 10" 2 


200 


a ~ 3.0 X 10" 


-4 


a ~ 10- 4 


300 


o ~ 3.5 x 10" 


-6 
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Figure 1. The Feynman diagram for the process e~ + d — > v e + \i~ \i~ + u through the X~ , W~ , 
and X gauge boson exchanges. This diagram shows that the doubly charged gauge boson X~ 
is produced through W~ and X exchanges between e and d and then decays into two muons. 
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Figure 2. The Feynman diagram for the process e~ + d — > v e + \i~ \i~ + u through the X~ gauge 
boson exchange. This shows that the dileptons are produced through W~ and X with the 
exchange of a neutrino. There are two graphs for this case. This is one of them. 
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Figure 3. The Feynman diagram for the process e~ + d — > v e + \i~ \i~ + u through the X~ gauge 
boson exchange. This shows that the dileptons are produced through W~ and X with the 
exchange of a neutrino. This is the cross term of the one in fig. 2. 
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Figure 4. The Feynman diagram for the process e~ + d — > v e + \i~ \i~ + u through the X gauge 
boson exchange. In this diagram, X is directly produced from the electron and decays into two 

muons. 



